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The power of second harmonic electromagnetic emission is calculated for the case when strong 
plasma turbulence is excited by a powerful electron beam in a magnetized plasma. It is shown 
that the simple analytical model of strong plasma turbulence with the assumption of a constant 
pump power is able to explain experimentally observed bursts of electromagnetic radiation as a 
consequence of separate collapse events. It is also found that the electromagnetic emission power 
calculated for three-wave interaction processes occurring in the long-wavelength part of turbulent 
spectrum is in order-of-magnitude agreement with experimental results. 
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Electromagnetic emission of a turbulent plasma at the 
fundamental plasma frequency lo p and its harmonics has 
been the subject of active theoretical and experimental 
research for several decades. This radiation carries infor- 
mation about properties of plasma turbulence, and reg- 
istration of plasma emission is one of the most efficient 
ways of studying physical processes occurring in space 
plasmas. That is why in most papers the problem of tur- 
bulent plasma emission is considered in the context of 
space phenomena such as type III solar radio bursts [T]- 
7 and emissions of planet's magnetospheres [EJEh Our 
interest to this problem is motivated by laboratory ex- 
periments |10j , in which electromagnetic radiation at the 
doubled plasma frequency is generated during injection of 
a high-current relativistic electron beam into the plasma 
confined in the GOL-3 multimirror trap. In contrast 
to previous beam-plasma experiments |11H14| . we study 
emission properties of rather hot (T =1-2 keV) plasmas 
in sufficiently strong magnetic fields (f2 = lu c /lo p ~ 1, 
where uj c is the electron cyclotron frequency). 

Several generation mechanisms of second harmonic 
plasma emission exist. In the context of type III radio 
bursts, either weakly turbulent coalescence of Langmuir 
waves £ + £ — > t [HI [16] , or generation of electromagnetic 
waves by collapsing caverns in strong plasma turbulence 
[51[T71[TH] are discussed. Moreover, such electromagnetic 
waves can be produced due to the Langmuir harmonic 
waves [19] scattering off density fluctuations. It is obvi- 
ous that the power of electromagnetic emission depends 
essentially on what nonlinear processes form the turbu- 
lent spectrum. In the theory of weak turbulence the im- 
portant role is played by the electrostatic Langmuir de- 
cay £ — > £' + s (s denotes ion-acoustic waves) and the 
reverse process £ + s — > £' . In an optically thick plasma, 
formation of a turbulent spectrum can be also affected 
by nonlinear processes £ —¥ t + s and £ + s — > t involv- 
ing electromagnetic waves with the plasma frequency u> p . 
In models of strong plasma turbulence [2"TH2"3"] wave en- 
ergy transfers through the spectrum due to Langmuir 
wave scattering off density fluctuations and collapse of 
localized Langmuir wave packets. Alternative models of 



strong plasma turbulence, in which wave collapse is su- 
pressed by either the direct conversion of Langmuir waves 
to damping modes [33] or radiative losses [35], are also 
discussed. 

Experiments on turbulent plasma heating at the GOL- 
3 multimirror trap show that intensity of electromagnetic 
emission demonstrates not only smooth variation in time, 
but also bursts with the duration of 2-10 ns, which we 
tend to associate with separate collapse events. Thus, the 
model of strong plasma turbulence, proposed in Ref. [21] 
and verified later by numerical simulations |26] . seems 
most appropriate for our experiments. In Ref. [3] this 
analytic model is used to estimate the power of electro- 
magnetic emission in the problem of type III radio bursts. 
In order to explain the results of laboratory beam-plasma 
experiments we will modify the model of Ref. [26| by tak- 
ing into account saturation of the pumping power due to 
beam trapping and generalize the calculation procedure 
of the second harmonic emission power of Ref. [3 to an 
arbitrary magnetic field. 

According to the model of strong plasma turbulence 
[2"Tl 121)] an isotropic turbulent spectrum of an unmag- 
netized plasma can be divided into three typical re- 
gions: source region, inertial range and dissipation re- 
gion. The source region occupies small wavenumbers 
k < fcjv/ — v"7 ( n T)/ r D (i"D is the Debye length) and 
consists of untrapped Langmuir waves, which are the 
products of beam-driven Langmuir waves scattering off 
long-wavelength density fluctuations. It is assumed that 
the spectral density of wave energy inside this part of the 
spectrum is independent on k. In the inertial range, mod- 
ulation instability results in trapping of Langmuir waves 
in local density wells and is followed by the wave collapse, 
which is responsible for the formation of the power-law 
spectrum. In the small-wavelength region various dis- 
sipation mechanisms come into force and spectral wave 
energy decreases rapidly with the increase of k. Thus, 
second harmonic plasma emission can be generated (i) 
due to Langmuir wave coalescence in the source region 
containing most of the wave energy and (ii) due to the 
collapse of trapped plasma oscillations which, despite the 
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low energy content, reach high energy densities at the late 
stage of collapse and result in radiation bursts. 

Let us calculate the level of the wave energy density 
W in a plasma turbulence that is pumped by an energy 
source with the power P, and estimate the typical du- 
ration of radiation bursts produced in separate collapse 
events. Energy balance between different parts of turbu- 
lent spectrum can be written in the form 
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The first equation describes the balance between pump- 
ing of beam-excited resonant waves with the energy den- 
sity Wr and dissipation produced by their scattering off 
long-wavelength density fluctuations. The rms level of 
these fluctuations is determined by the wave energy den- 
sity W concentrated in the source region: yj (6n 2 ) /n — 
aW/nT, where a = 0.7 is the numerical coefficient ob- 
tained in 2D simulations. The second equation in ([!]) 
shows that the power that comes to the source region 
from the pump is balanced by the power that leaves 
this region due to the wave collapse. Here, according 
to Ref. [55], the rate of collapse is determined not only 
by the rate of modulation instability, but also by the 
factor A(W) = 2\(W/nT) 1 / 2 accounting for finite time 
that takes the collapse to reach the self-similar regime 
(A ~ 0.7). 

The pump power is usually estimated as P = 2TWr, 
where T is the linear growth rate of the beam-plasma in- 
stability. In our experiments, however, the powerful elec- 
tron beam relaxes in the so-called trapping regime [27], 
for which the energy pump to the beam-excited coher- 
ent wave packets is saturated by the beam nonlinearity. 
Indeed, ID particle-in-cell simulations [2S] show that in 
our case the pump power does not depend on evolving 
parameters of the plasma turbulence. That is why in 
the balance equation this power can be considered as a 
given constant P = const. We can estimate this value 
from experimental data: P = f3nTo/ro = lOOkW/cm , 
assuming that it takes tq = 3 fis to heat the plasma with 
the density n — 2 • 10 14 cm -3 up to the electron tempera- 
ture To = 1 keV. The factor (3 = 6 takes into account that 
most of the dissipated wave energy goes to the formation 
of high-energy tails. 

For the resonant Wr and nonresonant W wave energies 
in this model we get 
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The characteristic duration of wave collapse can be esti- 
mated as 




For typical parameters of beam-plasma experiments at 
the GOL-3 multimirror trap, the wave energy density 
reaches the value W/nT = 0.01 for the electron tempera- 
ture T = 1 keV. The collapse duration at the same stage 
appears to be 3-4 ns, which is in a good agreement with 
the duration of radiation bursts observed experimentally. 
Thus, experimental data do not contradict the basic ideas 
of our model, and we can use estimates for the wave en- 
ergy density and the characteristic width of the energy- 
containing region of plasma turbulence to calculate the 
emission power. We also assume that these estimates are 
not extremely sensitive to the external magnetic field. 

Let us calculate the power of second harmonic plasma 
emission which is generated due to coalescence of un- 
trapped Langmuir waves t in the long- wavelength 

part of turbulent spectrum k < Um ■ By Langmuir waves 
in the cold magnetized plasma we mean plasma oscilla- 
tions pertaining to the upper-hybrid branch. Nonlinear 
interaction of such waves can be described in the frame- 
work of weak turbulence, but in contrast to the stan- 
dard theory we will take into account model damping of 
two-time correlation functions, which is used to describe 
the effect of finite life-time of Langmuir plasmons due to 
their scattering off density fluctuations. Let us represent 
the electric field in the following form 
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where uj^. and e£ denote eigenfrequencies and eigenvec- 
tors of linear plasma modes, E k — their slowly vary- 
ing amplitudes, and a indicates the branch to which 
they belong. In the general case, three-wave interaction 
a 1 + a" — » a is described by the equation 
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e(k, u>) is the dielectric tensor and y k is the Fourier 
transform of the second-order nonlinear electron current, 
which in the cold plasma limit takes the form 
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Figure 1. The spectral power dP/(dwd9) of x-mode emission for the isotropic turbulent spectrum (a) and for the spectrum 
including resonant waves (b). The spectral power dP / (dujd9) of o-mode emission for the isotropic turbulent spectrum (d) and 
for the spectrum including resonant waves (e). The spectral power dP/du) of x-mode (c) and o-mode (f) emission for the 
isotropic (dashed line) and anisotropic (solid line) spectrum. 
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In dimensionless units u p t, uj/lj p , xlo p /c, kc/ui pi 
eE k (uj p /c) 3 /(mcuj p ) for time, frequency, position, 
wavenumber and electric field amplitude of a plasma 
mode, respectively, three-wave interaction processes I + 
£ — > t between Langmuir and electromagnetic waves are 
described by the equation 



dEl 
dt 



1 



i i 



r<tu „ 

u k,k 1 ,k 2 x 



2(2n) 3 / 2 (dA i /duj) u 
e <(w£-<-<)* ( y(] £ _ kl _ k 2 )d 3 fcid 3 fc 2 . (8) 



According to Ref. [5] we assume that Langmuir 
waves, scattering off long-wavelength density fluctua- 
tions, change their phases stochastically with the char- 
acteristic frequency vjuj p = W /nT. In this case, the 
temporal correlation function of Langmuir electric fields 
can be written in the form 
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For the average energy of Langmuir turbulence we get 
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The spectral energy density of electromagnetic waves 
produced in spontaneous processes £ + £ — > t is governed 
by the equation 
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Thus, in the case of azimuthally symmetric turbulence, 
the spectral power of second harmonic electromagnetic 
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emission in units of nmc 2 is given by the integral 

o v ; 

where k(u>) is the solution of w = w[ and 8 is the polar 
angle of k. 

Let us compute the spectral power of electromagnetic 
emission for parameters typical for beam-plasma exper- 
iments in the GOL-3 multimirror trap. In the regime 
with the plasma density n — 2 ■ 10 14 cm -3 , the exter- 
nal magnetic field = 0.8 and the electron tempera- 
ture T = 1 — 2 kcV, the power of electromagnetic emis- 
sion was estimated experimentally as 0.1 -j- lkW/cm 3 . 
In our theoretical model, the isotropic part of long- 
wavelength plasma turbulence should occupy the spec- 
tral region k € (0.1; 2.45)w p /c and should contain the en- 
ergy W l /nT = 0.01 for the typical electron temperature 
T = 1 keV. Moreover, in the beam-excited turbulence 
there is an anisotropic population of resonant Langmuir 
waves containing about 5% of turbulence energy. This 
part of energy is concentrated in a rather small spectral 
region: k € (1.1; 1.3) and 9 £ (0;0.3). The computation 
results for the emission power of ordinary (o-mode) and 
extraordinary (x-mode) electromagnetic waves are pre- 
sented in Fig. [I] To analyze the contribution of resonant 
waves to plasma emission, we also present computations 
accounting for the isotropic part of turbulence only. 



The angular distribution of the emission power 
dP/dwdO shows that both x-mode [Fig. [l](a), (b)] and 
o-mode [Fig. [I] (d), (e)] are radiated predominantly in 
the transverse to the magnetic field direction. It is also 
shown in Fig. [Ijc) andjljf) that the total emission power 
integrated over angle and frequency reaches the value of 2 
kW/cm 3 and is dominated by the x-mode contribution. 
One can see that resonant waves do not substantially 
affect x-mode emission and result in the significant in- 
crease of the o-mode emission power. From the emission 
spectrum it is also seen that the main role in generation 
of electromagnetic waves is played by almost potential 
Langmuir waves with frequencies to > ui p . 

In conclusion, we calculate second harmonic electro- 
magnetic emission of a turbulent magnetized plasma 
driven by a powerful electron beam. We found that 
the simple analytical model of strong plasma turbulence 
with the assumption of constant pump power explains 
the results of laboratory beam-plasma experiments at 
the GOL-3 multimirror trap. We show that the typical 
duration of electromagnetic bursts observed in these ex- 
periments is in a good agreement with the theoretical es- 
timate of collapse duration. Our theory does also predict 
by order of magnitude experimental results for the total 
emission power and explains polarization of this emission. 

The author thanks Prof. A.V. Arzhannikov and Prof. 
LA. Kotelnikov for useful discussions. This work is 
supported by President grant NSh-5118.2012.2, grant 
11.G34.31.0033 of the Russian Federation Government 
and RFBR grants 11-02-00563-a, 11-01-00249-a. 



[1] V.L. Ginzburg, V.V. Zheleznyakov, Sov. Astron. 2, 235 

(1958) [15 

[2] D.A. Gurnett, R.R. Anderson, Science 194, 1159 (1976). 

[3] E.N. Kruchina, R.Z. Sagdeev, V.D. Shapiro, JETP Lett. [16 
32, 419 (1980). 

[4] R.J.M. Grognard, Sol. Phys. 81, 173 (1982). [17 
[5] M.V. Goldman, Sol. Phys. 89, 403 (1983). 
[6] PA. Robinson, I.H. Cairns, and A.J. Willes, Astrophys. [18 

J. 422, 870 (1994). 
[7] B. Li, PA. Robinson, I.H. Cairns, Phys. Plasmas 13, [19 

092902 (2006). 

[8] D.A. Gurnett, S.D. Shawhan, and R.R. Shaw, J. Geo- [20 

phys. Res. 88, 329 (1983). 
[9] I.H. Cairns and J.D. Menietti, J. Geophys. Res. 106, [21 

29515 (2001). 

[10] A.V. Arzhannikov, A.V. Burdakov, S.A. Kuznetsov, [22 
M.A. Makarov, K.I. Mekler, V.V. Postupaev, A.F. [23 
Rovenskikh, S.L. Sinitsky, V.F. Sklyarov, Fusion Sci. and [24 
Techno! 59 (IT), 74 (2011). 

[11] G. Benford, D. Tzach, K. Kato, D.F. Smith, Phys. Rev. [25 
Lett. 45, 1182 (1980). [26 

[12] P.Y. Cheung, A.Y. Wong, C.B. Darrow, and S.J. Qian, 
Phys. Rev. Lett. 48, 1348 (1982). 

[13] H.J. Hopman and G.C.A.M. Janssen, Phys. Rev. Lett. [27 
52, 1613 (1984). 

[14] A. B. Baranga, G. Benford, D. Tzach, K. Kato, Phys. [28 



Rev. Lett. 54, 1377 (1985). 

A. J. Willes, P.A. Robinson, and D.B. Melrose, Phys. 
Plasmas 3, 149 (1996). 

B. Li, A.J. Willes, P.A. Robinson, and I.H. Cairns, Phys. 
Plasmas 12, 012103 (2005). 

HP. Freund and K. Papadopoulos, Phys. Fluids 23, 732 
(1980); Phys. Fluids 23, 1546 (1980). 
K. Akimoto, H.L. Rowland, and K. Papadopoulos, Phys. 
Fluids 31, 2185 (1988). 

P.H. Yoon, R. Gaelzer, T. Umeda, Y. Omura, and H. 
Matsumoto, Phys. Plasmas 10, 364 (2003). 
V.N. Tsytovich, Theory of Turbulent Plasma (Consul- 
tants Bureau, New York, 1977). 

A.A. Galeev, R.Z. Sagdeev, V.D. Shapiro, and V.I. 

Shevchenko, Sov. Phys. JETP 46, 711 (1977). 

M.V. Goldman, Rev. Mod. Phys. 56, 709 (1984). 

P. A. Robinson, Rev. Mod. Phys. 69, 507 (1997). 

A.A. Galeev, R.Z. Sagdeev, V.D. Shapiro, V.I. 

Shevchenko, JETP Lett. 24, 21 (1976). 

W. Main and G. Benford, Phys. Fluids B 1, 2479 (1989). 

L.M. Degtyarev, R.Z. Sagdeev, G.I. Solov'ev, V.D. 

Shapiro, and V.I. Shevchenko, Sov. Phys. JETP 68, 975 

(1989). 

IV. Timofeev, K.V. Lotov, Phys. Plasmas 13, 062312 
(2006). 

IV. Timofeev, A.V. Terekhov, Phys. Plasmas 17, 083111 
(2010). 



